Origin of conductance quantization in disordered graphene ribbons 
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We present numerical studies of conduction in graphene nanoribbons with different types of dis- 
order. We find that even when defect scattering depresses the conductance to values two orders of 
magnitude lower than 2(? jh, equally spaced conductance plateaus occur at moderately low tem- 
peratures due to enhanced electron backscattering near subband edge energies if bulk vacancies are 
present in the ribbon. This work accounts quantitatively for the surprising conductance quantization 
observed by Lin et al. [Phys. Rev. B 78, 161409 (2008)] in ribbons with such low conductances. 
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Graphene nanoribbons (NR's) are attracting much 
experimentaliii^ii^ and theoretical^'^'^'^ '^°i^^'^^ interest 
due to their unique properties stemming from the Un- 
ear, massless Dirac-like spectrum of the underlying hon- 
eycomb lattice. In common with other quasi one- 
dimcnsional (ID) ballistic nanostructures, idea? NR's are 
expected to exhibit conductances quantized in integer 
multiples of the conductance quantum 2e^//i due to elec- 
tron transmission via subbands that arise from lateral 
confinement of electronic states in the NR. Recently Lin 
et ali- reported the first experimental observation of con- 
ductance quantization in NR's. However, surprisingly, 
the conductance steps that they observed were orders of 
magnitude smaller tha.n 2e^/h. They decreased in height 
with increasing NR length and were present only at mod- 
erately low temperatures lOK ^ T < 80K. Lin et aiji sug- 
gested that the quantized conductance steps that they 
observed may be due to different numbers of subbands 
in their NR's becoming populated with electrons as the 
back gate voltage in their system was varied. They at- 
tributed the low values of the quantized conductances 
(<< 2e'^/h) that they observed and their dependence on 
the length of the NR to low electron transmission prob- 
abilities through their device due primarily to scattering 
by defects. However, the quantum transport calculations 
reported to date^i^ i^"'^^ found conductance quantization 
to be destroyed by disorder even for NR's with much 
higher conductances (i.e., much less disorder) than those 
of the NR's studied by Lin et al} Still it should be noted 
that the theoretical work was for narrower NR's than 
those studied by Lin et al} and transport should be more 
sensitive to disorder in narrower NR's. Furthermore the 
possibility of quantized conductances in strongly disor- 
dered NR's has not been the subject of systematic the- 
oretical investigations. Thus the origin of the quantized 
conductances observed by Lin et alX has remained an 
open question. 

The purpose of this paper is to investigate theoreti- 
cally how different scattering mechanisms affect electron 
transport in wide disordered NR's such as those of Lin 
et air and to clarify under which conditions quantized 
conductances much smaller than 2e^/h can occur in such 
systems. We find that electron scattering by C atom va- 
cancies at moderately low temperatures (in the presence 



of edge disorder) can account quantitatively for the quan- 
tized conductances observed by Lin et al} The underly- 
ing mechanism that we identify is modulation of the NR 
conductance by enhanced electron back-scattering by va- 
cancies whenever a subband edge crosses the Fermi level. 

We concentrate on three disorder types, namely bulk 
vacancies, edge imperfections and long-range potentials 
due to charged impurities. Other disorder types such as 
weak short-range potentials due to neutral impurities,— 
and lattice distortionsAiI, may be present, but the three 
disorder cases to be discussed here have the strongest im- 
pact on transport through NR's and thus are more rel- 
evant to the strong conductance suppression reported in 
Ref. ill. Among the three, only bulk vacancies scale the 
heights of different conductance steps uniformly. Thus 
their presence appears crucial for the observation of con- 
ductance quantization in NR's with strong disorder. 

We describe NR's by the standard tight-binding Hamil- 
tonian on a honeycomb lattice, 

where et is the on-site energy and tij = t = 2.7 eV is the 
matrix element between nearest-neighbor atoms. This 
Hamiltonian is known to describe the it band dispersion 
of graphene well at low energies. Spin and electron in- 
teraction effects are outside of the scope of our study. 
Bulk vacancies and edge disorder are introduced by ran- 
domly removing carbon atoms and setting appropriate 
hopping elements tij to zero. It is assumed that atoms 
at the edges are always attached to two other carbon 
atoms and passivated by a neutral chemical ligand, such 
as hydrogen. The bulk and edge disorder are charac- 
terized by the probability of the carbon atoms being re- 
moved, and p'^, respectively, p^ is normalized relative 
to the whole sample, while p'^ is defined relative to an 
edge only. The long-range potential due to charged im- 
purities is approximated by a Gaussian form^'^^ of range 
d: ei — J2ro K)exp(— \ri — ro|^/d^), where both the am- 
plitude Vq and coordinate tq are generated randomly. 

In the linear response regime the conductance of the 
NR is given by the Landauer formula)^ 

^ 2e^ r ^.^^.dfME'-E) 
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FIG. 1: (color online) Conductance (a)-(c), average conductance (d)-(f) and conductance fluctuations (g)-(i) as a function of 
energy for the graphene ribbon of width W = 30 nm and different lengths L — 10. ..1000 nm. Left panel corresponds to bulk 
vacancies, middle panel is for edge disorder while right one shows the effect of long-range potential. Parameters of disorder 
with representative illustrations are given in the plots (a)-(c). The dotted lines in (a)-(f) show the conductance quantization 
for the ideal ribbon. The gray filled areas in (d)-(f) denote DOS for the ideal ribbon. The dashed lines in (g)-(i) mark the 
universal value of the conductance fluctuations for quasi-lD systemsiii Temperature T = 0. t = 2.7 eV; a = 0.142 nm. 



T{E) is the total transmission coefficient and fpoiE) 
is the Fermi-Dirac function. T{E) is calculated by the 
recursive Green's function method, see Ref. 12 for de- 
tails. Fluctuations of the conductance are defined by 
rms{G) = {{G'^) — {G)^Y^'^^ where () denotes averaging 
over an ensemble of samples with different realizations of 
disorder. For the results presented below, averaging was 
carried out over ten realization for each disorder type. 

To investigate the transport properties of disordered 
NR's we choose geometries similar to ones studied 
experimentallyii The disorder is assumed to exist in a 
finite ribbon of width W and length L. This ribbon is 
attached at its two ends to semiinfinite leads represented 
by ideal NR's of width W . The edge (host) configuration 
is taken as armchair in the following. Representative dis- 
order geometries are shown in the insets in Figs. [Ua)-(c). 

Fig. [1] shows the effect of different disorder types on 
conduction in NR's. For each disorder type we keep the 
defect concentration and strength and the ribbon width 
fixed {W — 30 nm) and vary its length L. As L increases, 
and the number of scattering centers grows, the conduc- 



tance decays and quantization steps are destroyed. 

For bulk vacancy disorder, even a small concentration 
of the defects affects the conductance strongly; see Figs. 
[Ua),(d). Apart from reduced conductances, the disor- 
der results in sample-specific conductance fluctuations, 
Figs. [Ha),(g), whose amplitude is of order e^//i, inde- 
pendent of energy or NR length. This is a quantum in- 
terference effect similar to the universal conductance fluc- 
tuations (UCF's) of mesoscopic metalsJ^ The particular 
value of the conductance depends sensitively on the elec- 
tron energy, ribbon length and locations of the vacancies. 
Since the vacancies are distributed over the whole sam- 
ple, intra-subband scattering predominates. Thus the 
conductance in Figs. [Ua),(d) (coarse grained in energy 
to smooth out UCF's) scales uniformly with NR length 
L, i.e., in a similar way for all subbands. This resembles 
bulk island scattering in conventional quantum wiresJ^ 

By contrast, for edge disorder in Figs. [T]Ib),(e) the 
conductance scales non-uniformly: Defects at the bound- 
aries scatter electrons equally into all subbands result- 
ing in stronger suppression of the conductance at higher 
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energies E where more subbands are available; see e.g., 
L = 1000 nm ribbon in Figs. [IIb),(e). 

Potential inhomogeneities due to charged impurities 
lead to the appearance of electron and hole puddles in 
NR's.^* Scattering by the potential inhomogeneities re- 
sults in subband mixing that smears conductance steps, 
Figs. [Hc),(f). As the subband number increases, in- 
tervalley scattering becomes more effective with stronger 
backscattering of higher subband states in long ribbons. 
The first subband, however, is not affected by the long- 
range potential because of internal phase structures of 
its wave function that make the scattering amplitude 
vanishiii The conductance fluctuations are roughly twice 
as strong as for bulk vacancy and edge disorder, Figs. 
[Ijg)-(i). This may be due to weaker inter- valley scat- 
tering for which particles at K and K' Dirac points 
contribute independently to the UCFs.— The fluctua- 
tion amplitudes agree reasonably well with the value for 
UCF's in quasi-lD systemgi^ 0.729^. 

A prominent effect of all disorder types is the forma- 
tion of a conductance dip when the Fermi level crosses 
a subband edge. This is most obvious in the averaged 
conductance (G), Figs. [Tlld)-(f). The origin is the strong 
intersubband scattering caused by defects, where an elec- 
tron in a state \nk) scatters into another state \n'k'). It 
can be understood physically by considering the Fermi 
Golden rule expression for the scattering time r— : 
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j-J2\{nk\H'\n'k')fp„, {E) . 



(3) 



Here H' is the perturbation due to defects and pn' [E) 
the density of states of the n'-th subband. Assuming that 
\{nk \H'\ n'k')\^ is independent of the band index n', the 
scattering rate 1 /r is seen to be proportional to the total 
density of states of the ribbon, p{E) = pn'{E). For 
a perfect ribbon the dispersion relation can be approxi- 
mated by a parabolic function if k is small and |n| > 1.— 
Therefore, p{E) diverges at subband thresholds En' as 
{E-En'Y/^. This agrees with the numerically calculated 
density of states for the tight-binding Hamiltonian ([1]), 
see the gray areas in Figs. [Hd)-(f). Thus, the scattering 
time r is strongly reduced when the Fermi energy ap- 
proaches a subband threshold En' and the transmission 
of electrons in the n-th subband is strongly suppressed 
due the scattering into the other n' subbands. As a re- 
sult, the conductance shows dips at the subband edges. 

As the temperature increases, the conductance fluctu- 
ations are smeared out and the dips in the conductance 
associated with enhanced electron back scattering when 
the Fermi level crosses subband edges become clearly visi- 
ble for temperatures T not greatly exceeding the subband 
energy separation, AirksT w AE = En+i — En- For 
graphene ribbons 30 nm wide AE w 0.02/; = 54 meV, see 
Fig. [H that corresponds to T « = 50 K. Above 

this temperature the conductance dips become gradu- 
ally smeared, but well below it the conductance may be 
dominated by UCF's of the disordered ribbon. This es- 
timate is in good agreement with calculations presented 
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FIG. 2: (color online) Conductance through disordered rib- 
bons as a function of the Fermi energy for temperatures 
T = 0, 80, 300 K. The ribbons have width W = 30 nm and 
length L = 1000 nm; parameters of disorder are listed in Fig. 
[T] Thin solid vertical lines correspond to energies when num- 
ber of subbands changes by one. For the sake of clarity the 
curves for long-range potential are shifted upward by 2e^/h. 



in Fig. [2l where ribbons with different disorder types are 
subjected to T = 0,80,300 K: The conductance dips at 
subband edges manifest as the smooth conductance os- 
cillations that are clearly visible for T = 80 K. They are 
very regular and are superimposed on a smoothly rising 
background for the case of the bulk vacancies, but ap- 
pear very distorted when the edge disorder or long-range 
potential introduced, except for the first two or three os- 
cillations for the case of the edge disorder. Since the 
effects of long-range disorder are similar to those of edge 
disorder (see Figs. [ire),(f)) we shall not consider the 
long-range disorder further here. 

The conductances of the NR's measured by Lin et alX 
were 65-260 times smaller than the conductance quan- 
tum and also much smaller than the conductances of the 
model systems studied above in Figs. [1] and [21 How- 
ever, the ideas developed above apply equally well to the 
lower conductance regime in which the experiments were 
carried out and are able to account quantitatively for 
the conductance quantization that Lin et alX observed. 
We demonstrate this next by presenting simulations for 
NR's with the same sizes as in the experiments— and with 
defect concentrations chosen to yield low conductances 
similar to those measured by Lin et alii 

Fig. [3] shows the calculated conductances (G) of the 
disordered ribbons along with the experimental data 
from Ref. [l|. The features in the theoretical plots that 
match the experimental conductance plateaus are the 
conductance dips that are due to enhanced electron back 
scattering at the energies of the subband edges of the 
nanoribbon that we have already discussed in connec- 
tion with Figs. [Hand [21 The agreement between theory 
and experiment is remarkable especially for the heights of 
the conductance plateaus. From the theoretical point of 
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FIG. 3: (color online) Comparison of theoretical and experi- 
mental data for the ribbons oi W — 30 nm and L = 850, 1700 
nm. Experimental data is adopted from Ref. [l|. Theoretical 
calculations are performed for two disorder cases: (a) - only 
bulk vacancies p'' — 4 x 10~^ (solid red curve), and combi- 
nation of edge p'^ = 2 and bulk p'' = 10"'* disorder (dotted 
green curve); (b) - p*" = 8 x 10"*, and combination of p" = 3.5 
and — 2 x 10^*. The gate voltage Vg is scaled to produce 
better fit; note that E ~ ^ discussed in Ref[l9l. 

view, there are several detailed scenarios that might re- 
sult in this behavior. Unfortunately, it is not possible to 
rule out any of them because the experiment gives no in- 
formation regarding which disorder is type actually real- 
ized. Therefore, we propose that the dominant scattering 
mechanism might be either due to bulk vacancies alone or 
a combination of rough edges with a lower concentration 
of bulk vacancies. The presence of the latter is crucial 
because they equalize the differences between the con- 
ductances of the different plateaus making them equidis- 
tant. In particular, we found that p*" = 4 x 10""' bulk 
vacancies are enough to reduce the conductances of the 
quantized plateaus by a factor of 65 relative to the con- 



ductance quantum in accord with the experimeniji, 
see solid red line in Fig. ^a,). This means that one 
in 2500 carbon atoms is removed, which seems plausi- 
ble. The other scenario consists of distorted edges with 
two rows of carbon atoms removed on average along the 
boundaries and also one in 10000 bulk carbons removed, 
p'^ = 2 and p"^ — 10"^, see dashed green line in Fig. 
[3fa). For the longer L — 1700 nm ribbon the height of 
conductance steps drops to a factor 260 lower than the 
conductance quantum. Fig. [SJb). This implies defect 
concentrations twice those of the shorter L = 850 nm ex- 
perimental ribbon. The lower temperature in Fig. [Sj^b) 
results in stronger conductance fluctuations than in Fig. 
^a); the fourth plateau being not discernible in the ex- 
perimental data'^ in Fig. El^b), may also be due in part to 
a particular disorder configuration. However, all visible 
conductance plateaus are due to subband formation asso- 
ciated with particle motion quantized in the transverse 
direction. At much lower temperatures in our simula- 
tions these conductance plateaus are not discernible due 
to UCF's and they also disappear completely at room 
temperature, behavior similar to that in Fig. [21 and com- 
pletely consistent with the data of Lin et alX. 

In conclusion, our quantum transport calculations 
have shown that equally spaced quantized conductance 
plateaus should be observable in disordered graphene 
nanoribbons even for conductance values much smaller 
than the conductance quantum 2e^ /h at temperatures 
comparable to subband energy spacings. The plateaus 
are due to enhanced electron back scattering by defects 
at energies near subband edges. Deviations from equal 
spacing of the conductance plateaus can occur depend- 
ing on the defect configurations in particular experimen- 
tal samples. These findings provide a microscopic ex- 
planation of the conductance quantization of graphene 
nanoribbons observed by Lin et alX and suggest that the 
observed conductance quantizationi can be regarded as 
a signature of subband formation. 
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